Mean p t fluctuations from 2- and 4-particle correlations 
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Recent results on event-by-event mean transverse momentum, (pt), fluctuations in ultra- 
relativistic heavy ion collisions are briefly reviewed. We conclude that the observed fluctuations 
are in a rough agreement with that expected for the independent superposition of nucleon-nucleon 
collisions. We further discuss the possibility of the use of the forth order cumulants of the particle 
transverse momentum distribution in order to access the fluctuations related to collective phenom- 
ena. 
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One of the most interesting questions in the field of 
the ultra-relativistic heavy ion collisions is the question 
of the hadronization of the system. Is the system ther- 
malized/equilibrated? Does the system evolution include 
the phase transition? Event-by-event fluctuations, e.g. 
of the mean transverse momentum, are considered to be 
one of the important tools to answer these questions jlj . 
The fluctuations depend on the nature of the phase tran- 
sition. A first order phase transition may lead to large 
fluctuations in energy density due to formation of QGP 
droplets^], || . Second order phase transitions may lead to 
divergence in the specific heat; it would also increase the 
fluctuations in energy density due to long range correla- 
tions in the system |I| . One could observe them as fluc- 
tuations in mean transverse momentum if matter freezes 
out at the critical temperature T c j|, ||, Q, Q|. 

The centrality dependence of the fluctuations is an im- 
portant observable. If the fluctuations are due to the 
particle production via some kind of clusters (e.g., res- 
onances, strings, (mini)jets, independent TViV-collisions, 
etc.) and the relative production of such clusters do not 
change with centrality, the magnitude of the dynamical 
part in fluctuations should be inversely proportional to 
the number of the clusters, therefore to the particle mul- 
tiplicity. New physics should appear as a deviation from 
such a dependence. There could be two reasons for the 
change in the centrality dependence. First, such collec- 
tive phenomena as phase transition affect many particles 
in the system (unlike as in the scenario of particle pro- 
duction via a few particle clusters), and, therefore have 
different multiplicity (centrality) dependence. Second, 
new phenomena are expected to happen at some critical 
particle density, which in turn depends on centrality. 

Consider fluctuations in the inverse slope of mt = 
\J m 2 + p\ distribution (effective temperature fluctua- 
tions). Such fluctuations, for example, could be due to 
event-by-event fluctuations in radial expansion velocity. 
Then, depending on the slope, on average, the transverse 
momenta of all particles would be lower or higher com- 
pared to the average over all events. It results in (posi- 
tive) correlations between transverse momenta of differ- 
ent particles ]§, M. One can quantify such correlations 



by the particle transverse momenta covariance 

K 2 ,p t = COv(p t ,i,Pt,j) = (SptiSPtjhitj = 2 { Pt ),dynam> i 1 ) 

where 5pa — pt.i — pi , with pi being the inclusive mean 
transverse momentum. It can be also useful to consider 
dimensionless quantities like rf Pthdyna J <$ t , 
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from the fact that this quantity equals to the difference 
between the actual variance of the event mean p t dis- 
tribution and the expected width due to the statistical 
fluctuations Jf9| in (p t ). The latter is due to the finite 
event multiplicity. 
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and M is the multiplicity. 

Event-by-event dynamical fluctuations have also been 
analyzed by several experiments using the so called 
$ pt | 10 measure (the approximate relation to a 2 (pt) dynam 
is taken from ||): 



$p t = V(((Pt) - Mp t ) 2 )/(M) - (T pu inclv 
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and very close to it the difference factor Aer Pt [ pd| : 

^^Pt — ( \J (-^-0 &Pt ^inclusive @ pt .inclusive) ~ *^r, 



Using the above relations to compute K2,p t whenever 
$ pt being reported, we arrive to the conclusion that 
in central collisions of heavy nuclei such as gold and 
lead, relative fluctuations in mean transverse momen- 
tum, ^(p^^dynam/Pt is of the order of 1 - 1.5%. For the 
6% most central collisions STAR reports |ll|] the prelim- 
inary results of <r {pthdynam /pi = 1.2 ± 0.2. The PHENIX 



2 



measurements [|12| have larger uncertainties, but if av- 
eraged over different centralities (assuming <& pt does not 
change with centrality), this measurements yield for the 
central Au+Au collisions the value of cr/p t \ ^y nam /pl = 
1.4 ± 0.9. In Pb+Pb and Pb+Au collisions at CERN 
SPS (Jsnn = 17 GeV) 3>p t have been measured by 
NA49 p| and NA45/CERES jl§]. The published NA49 
result |13| & Pt = 0.6 ± 1 MeV/c was obtained for the 
rapidity region 4.0 < y < 5.5 and 5% most central col- 
lisions. The observed fluctuations are extremely small, 
but it would be incorrect to compare these numbers with 
RHIC measurements, since NA49 results were obtained 
in the forward rapidity region. The CERES collaboration 
at the SPS has measured the fluctuations in the central 
rapidity region. They report <& Pt = 7.8 ± 0.9 MeV/c 
in central Pb+Au collisions [Q. As seen from Eq. 
$ Pt is directly proportional to the number of recon- 
structed tracks used for its calculation (subject to ac- 
ceptance cuts and tracking efficiency), which compli- 
cates the comparison of the results from different experi- 
ments. For a rough comparison one can take into account 
the CERES multiplicity ~130, ((p t )> «420 GeV/c and 
cr pt , inclusive « 0.270 GeV/c. Then $ Pt w 8 MeV corre- 
sponds to a ipA .dynam /¥t ~ 1.4%. The centrality depen- 
dence of (pt) fluctuations, whenever studied, is consistent 
with particle production via clusters picture within about 
30%. 

The dynamical part of mean transverse momentum 
fluctuations has been measured at the ISR JT5[ in pp 
collisions. This was done by analyzing the multiplic- 
ity dependence of <J( Pt ) under the assumption that in 
pp collisions the dynamical part in (p t ) fluctuations 
does not depend on multiplicity. It was observed that 
0-<p t >, dynam/ '((Pt)) ~ 12 %- Rescaling of this quantity 
with (the square root of) the ratio of multiplicity den- 
sities in pp and Au + Au collisions yields the fluctuations 

in Au + Au: CF ( Pt ) ^dynam 

/((pt)} w 0.8%, about 50% less 
than that observed in AA collisions. Rescaling the quan- 
tity with the number of participants gives even better 
agreement. 

The non-zero value of H2.p t could be due to different 
reasons: including such as resonance decays or jet pro- 
duction, which are not "interesting" from the point of 
view of the collective phenomena under search in nu- 
clear collisions. In principle, selecting specific pairs for 
the averaging in Eq. [I] one can suppress or enhance the 
contribution to (j>t) fluctuations of different origin. For 
example, one can correlate: (1) particles in two disjoint 
pseudo-rapidity regions. The 'gap' between the two re- 
gions eliminates effects such as quantum statistics (Bosc- 
Einstcin and Fcrmi-Dirac) correlations or Coulomb final 
state interactions. (2) Positive particles with negative 
particles, which is expected to enhance the contribution 
from resonances. All those tricks do not solve the prob- 



lem completely. Below we discuss the possibility of the 
use of 4-particle correlations (cumulants) in order to eval- 
uate genuine collective phenomena contribution to mean 
p t fluctuations: 

K 4,p t = (SptiSptjSptkSptm) - ZiSptiSptj) 2 (5) 

with all i, j, k, and m being different. The relative con- 
tribution of the collective phenomena and particle pro- 
duction via clusters to the cumulant K4 tPt is enhanced by 
a factor of ~ M tota i / 'M c i uste r compared to the relative 
contribution to H2.p t ■ The notation M c i uster is used here 
for the average cluster decay multiplicity, and M to tai is 
the total event multiplicity. 

In jl6) cumulants have been proposed for the study 
of multiplicity fluctuations. Recently, four particle az- 
imuthal correlations (cumulants) jlTj have been success- 
fully used by the STAR Collaboration in the analysis of 
the azimuthal correlations for the measurements of el- 
liptic flow |l8| |. With modern statistics the technique 
works well even for relatively small signals (two particle 
correlations of the order of 10~ 4 and smaller. This is 
just the range of the correlation magnitude needed for 
the 4-particle correlations mean p t analysis (recall that 
K 2.p t /{pt) ~ (1 — 1.5%) 2 . In many particle cumulant 
mean p t analysis one can use the generating functions 
similar to the ones used in (It], |l8) . 

The contribution to Ki pt from the inverse slope fluc- 
tuations to the first order is proportional to the cor- 
responding cumulant of the inverse slope distribution. 
For the case of of an equal mixture of event with two 
different slopes T± and T2, the corresponding cumulant 
k 4 ,t = -(AT) 4 /8, (AT = Ti - T 2 ). Note that for some 
specific inverse slope distributions the cumulant could be 
small (for a rectangular distribution with width AT it 
is K4 : t = — (AT) 4 /120, and for a Gaussian distribution 
it is zero). The very small values of the cumulants may 
complicate the analysis. The use of additional weights in 
Eq. @ could somewhat help in this case. 

In summary, we observe that while clear correlations 
in particle transverse momenta has been measured in 
heavy ion collisions at SPS and RHIC, no evidence that 
the correlations are due to collective phenomena have 
been found so far. The centrality dependence of the ob- 
served correlations, as well as comparison to pp collisions 
at similar energies, is roughly consistent with the picture 
of particle production via clusters. In order to disen- 
tangle the contribution to the mean p t fluctuations com- 
ing from collective phenomena we propose to use many- 
particle correlations. 
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